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Abstract. By using dark nebulae as opaque outer screens, the 
luminosity function of white dwarfs and M dwarfs has been 
studied. High-extinction areas towards the Orion A, Serpens 
and p Oph cloud complexes were observed, covering a volume 
corresponding to 464 pc'' in the solar neighbourhood, complete 
to A/v ^ 16.5. Foreground stars were selected by VRI pho- 
tometry and photometric parallaxes. 

The resulting foreground sample consists of 21 M dwarfs 
and 7 white dwarfs. The derived M-dwarf luminosity function 
is consistent with previous studies, showing no substantial up- 
turn beyond Afy = 16. The 7 white dwarfs imply a local 
density of 0.013 (±O.OO5)M0 pc^^ ^5% (he dynami- 
cal mass in the solar neighbourhood) supporting other recent 
results but obtained with a completely different method. 

For the clouds, foreground-star data were used to derive 
independent cloud distances, while the background stars and a 
simple model of the Milky Way gave reliable cloud extinctions. 

Key words: Stars: late-type - low-mass, brown dwarfs - lumi- 
nosity function, mass function - white dwarfs - Galaxy: solar 
neighbourhood - stellar content 



1. Introduction 

The dark matter affects astronomy on all distance scales, single 
galaxies, groups of galaxies and the Universe as a whole. The 
dynamical mass in the solar vicinity has recently been deter- 
mine d from Hipparcos data to 0.076 ± 0.015 Mq pc^"^ (Creze 
et al. 1997), leaving little space for local dark matter, since the 
visible matter in form of stellar-like objects and the interstellar 
medium comprise ^ 0.08 Mq pc~'^. 

The stellar contribution to the local mass density derived 
by Jahreiss & Wielen (1997) is 0.039 Mq pc~'^, which agrees 
within errors w ith mo st other recent figures (e.g. 0.0 42 ± .008 
in Mera et al. ( |1996| ), 0.05 ± 0.01 in Kroupa et al. ( |1993| ) and 
0.034 in Gould et al. ( 1996 )). Similar amounts of matter have 
been suggested to reside in the interstellar medium, ~ 0.04 



M q pc~ -^ (Haywood et al. |1997| ; Creze et al. |1997| ; Bahcall et 
al. 1992 and references therein). However, present estimates of 
the local surface density of the interstellar medium range from 
6 to 13 Mqpc~^, indicating a substantial uncertainty also in 
the space density. 

The photometric surveys by Gould et al. ( 199^ and Kirk- 
patrick et al. (1994) give 0.01 1 Mq pc^^ for M dwarfs in the 
mass interval 0. 1-0.5 Mq. For the nearby stars the correspond- 
ing numbers are 0.013 and 0.015 Mqpc~^ respectively for 
Reid & Gizis (|1997|) (hereafter RG97) and Jahreiss & Wie- 



len (1997). The discrepancy between the photometric and the 
nearby samples is significant and has been suggested to arise 
from unresolved binarity (Kroupa 1995) and structure in the 
Afy vs F — / relation that has not been taken into account 
(RG97). For the white dwarf (WD) mass density, the 0.003 
WDspc-3 (ro 0.002 M0PC-3) given in Liebert et al. ( [198^ 
(hereafter LDM88), is less than 3 % of the local dynamical 
mass. There are, however, recent investigations that p oint to 
a much higher mass hidden in WDs. Ruiz & Takamiya ( 1995| ) 
(hereafter RT95) derived, from a large proper motion survey, 
a factor three higher space density t han in LDM88. A similar 
value was obtained by Oswalt et al. ( 1996 ) from WDs in bina- 
ries. This revision of the WD space density (supported also by 
this paper) indicates that as much as 15 % of the local dynami- 
cal mass can be hidden in WDs. 

Thus, although the bulk of the local mass seems to be iden- 
tified, the role of the late M dwarfs and WDs in this context is 
still not clear and remains a question to discuss. 

In this paper the subject is addressed in a survey of opaque 
nebulae, obtaining volume limited and dynamically unbiassed 
samples of foreground M dwarfs an d WD s. This method was 
first proposed by He rbst & Dickm an ( 1983 ) and h as prev iously 
been used by Jarrett ( |l992|Jl995| ) and Jarrett et al. ( |1994|) (h ere- 
after JDH94), whose results will be discussed in Sect. 8.1. The 
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* Based on observations collected at the European Southern Ob- 
servatory (ESO), La Silla, Chile and at the Nordic Optical Telescope 
(NOT), La Palma 



present survey covers a volume corresponding to 464 pc^'^ in 
the solar neighbourhood and is complete to Afy ^ 16.5. The 
completeness limit is 2 magnitudes deeper than the previous 
surveys by JaiTett and the volume twice as large. 

In § 2 the selection of the target clouds is described. § 3 out- 
lines the observations and reduction procedures. The extinc- 
tions of the clouds are derived in §4. The selection of fore- 
ground stars, including colour-colour analysis and photomet- 
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ric parallaxes is decribed in § 5. Constraints on our foreground 
sample obtained from other sources are given in § 6. The M- 
dwarf luminosity function (LF) and WD space densities, in- 
cluding corrections for multiplicity and galactic density gradi- 
ents are derived in § 7. In § 8 the results of this paper are com- 
pared to other recent findings, with special emphasis on the 
dark cloud survey described in JDH94. The conclusions and 
some ideas for the future are given in § 9. 

2. Selecting target clouds 

A suitable target cloud should fulfill certain criteria. Its dis- 
tance has to be known in order to match integration times to 
the faintest stars at the cloud's edge. It has to be sufficiently 
opaque (Ay <; 5 mag) to screen out the vast majority of back- 
ground stars and galaxies. The high extinction areas must be 
large enough to provide a reasonable number of foreground 
stars. It is also desirable to avoid the immediate neighbourhood 
of star-forming regions, since young overluminous stars within, 
or just in front of, the cloud may be mistaken as nearby older 
field M dwarfs. 

The final selection includes the following cloud complexes, 
Orion A (southern part of LI 641, 480 pc, / = 213°, 5= -19°), 
Serpens 4° south from the Serpens cloud core, 200-250 pc, 
/ = 27°, b = 3°) and two separate areas in p Oph (L1688 and 
L1689, 160 pc, I = 354°, b = 16°). Cloud distanc es for the 
original selection were taken from Hilton & Lahulla (1995). 



2.7. Cloud distances 

The photometric distances and the total number of M dwarfs in 
front of Orion A indicate a distance of 400 rather than 480 pc. 
This result is consistent with the 390 pc in Anthony-Twarog 
(1982), derived from a subsample in the Warren & Hesser 
(1978) catalogue of stars associated with the Orion complex. 
In order to confirm this finding, the Hipparcos catalogue (ESA 



1997) was searched, resulting in parallaxes with standard er- 



rors less than 1 mas for 92 of the Warren stars, showing a well- 
defined peak at ~ 400 pc. 



For Serpens, Straizys et al. (1996) gave a mean distance of 
259 pc for 18 stars associated to the cloud. Their total sample 
consists of 99 stars, foreground, embedded and background. 
The low-extinction part of their sample indicates an extinction 
wall at ^ 255 pc. Of the 99 stars, 13 were measured by Hippar- 
cos, revealing a few cases of misclassified luminosity classes, 
but basically confirming a distance of ~ 260 pc. 

Eleven of the stars that were assigned a high probability to 
be associated with the p Oph complex by Elias (1978) were 
identified in the Hipparcos catalogue and show a sharp peak at 
150 pc. 

The distances finally adopted in this paper are 400, 255 and 
150 pc respectively to the Orion A, Serpens and p Oph com- 
plexes. Note that these clouds are extended irregular objects far 
from being sheet-like or spherical. The distances above were 
derived from observations of optically visible stars, and there- 
fore mainly probe the front edge of the cloud, which though is 



sufficient for the purpose of this paper The errors in the cloud 
distances are expected to be on the order of 10 %, giving a cor- 
responding volume uncertainty of ~ 30%. 

3. Observations and reductions 

The main body of data consists of VRI photometry and was 
collected at the ESO NTT 3.55 m telescope at La Silla, Chile. 
The Serpens and p Oph areas were observed in Jun28-Jul01, 
1995 and cover 536 and 816 arcmin^ respectively. The Orion 
A field, observed in Jan24-26 1996, covers 158 arcmin^. The 
same instrumental setup was used in both runs, EMMI RED 
with CCD #36, TEK2048, each pixel corresponding to 0.27". 
The unvignetted field of view was 8.6' x 9.15'. Bias and flat- 
field corrections were done in a standard way withlRAFQ. 

The transformations to the Johnson V and Kron-Cousins 
RI sy stems were based on standard stars selected from Landolt 
(1992). Since the observations were aimed at M dwarfs in the 
first place, care was taken to include the reddest dwarf stars in 
this list, G45-20 (V ~ I = 4.00), G156-31 (V - I = 3.68), 
G12-43 (V - I = 3.48), GlO-50 (V ~ I ^ 2.97) and G44- 
40 {V — I = 2.79). The transformation to the standard system 
includes corrections for extinction, zero-point offsets and linear 
colour terms. The rms of the standard star residuals across the 
whole fitted colour range (0 < — / < 4) is 0.02 mag in VR 
and 0.03 mag in /. 



3.1. Complementary NOT observations 

For Orion A, part of the VRI photometric candidates were 
measured in J using ARNICA, a NICMOS3 near-IR array at 
the Nordic Optical Telescope (NOT), La Palma, in Sep 1996. 
The J calibrations were made via standard stars provided by 
the ARNICA team (Hunt et al. |1995|). The Hunt transformation 



was adopted: Jcit — '/arnica+ zero-point offset. The final 
1 cr error in the calibration of the J magnitudes is 0.07 mag. 

Low-resolution spectra of the two brightest WD candidates 
were obtained at the NOT in Nov 1997, using the ALFOSC, a 
low-resolution spectrograph/focal reducer combination instru- 
ment. Standard reductions were carried out with IRAF. It was 
not possible to apply a meaningful flux calibration, since the 
targets were always at airmasses between 2 and 3.5, leading to 
a substantial loss of light in the blue end of the spectra. 

Since this survey was designed for detecting M dwarfs, the 
original / exposure times were set rather short. The discov- 
ered WDs are bluer than the M dwarfs, resulting in uncertain 
/ magnitudes in the original data. The WD / photometry was 
improved by ALFOSC observations at the NOT in Nov 1997. 
Again standard methods were used in the reductions. 

All instrumental magnitudes were extra cted b y an empiri- 
cal growth-curve method outlined in Festin (1997). 



IRAF (Image Reduction and Analysis Facility) is distributed by 
National Optical Astronomy Observatories (NOAO), which is oper- 
ated by the Association of Universities for Research in Astronomy, 
Inc., under contract with the National Science Foundation. 
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Fig. 1. My at the peak in the through-shining background 
dwarf distribution as a function of the cloud's visual extinction. 
The solid line is a polynomial least-square fit to the calculated 
points (asterisks). The rings mark the three reference fields at 
the extinction Ayi as derived in Fig. ^ 
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Fig. 2. The two-step extinction estimation. The dashed line is 
the main sequence. The solid line is the GOV reddening curve 



in Jarrett (1992). Each dot is a 0.5 mag step in Ay. The arrowed 
solid lines show the dereddening steps. The rings demonstrate 
the location of a fictious set of stars 



The My completeness limit of the survey was set by the V 
data, as the magnitude at which the mean of the internal mag- 
nitude errors in a field equals 0. 12 mag. This limit was justified 
by a star-count analysis in a previous study of the Pleiades (Fes- 
19981 ). Table [l| es completeness hmits and positions for 



tin 



all the observed fields. 



4. Star counts and cloud extinctions 

The extinction offsets relative to a reference field were derived 



for each field by the method of star counts (Dickman 1978). 
Since there are more background stars present in the / fields 
than in V or i?, the / band was used for the star counts. The 
derived extinction in / is transformed to V by the relation in 
Winkler (|l997|), Ay ^ 1.69Ai. 
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Fig. 3. The extinction in the reference fields. The main se- 
quence and reddening curves are as in Fig. |[ The rings are 
real stars and the horizontal bars show steps 1 & 2 of Fig. |2| 



In order to estimate absolute extinctions in the reference 
fields, a two-step method was used. A first estimate is ob- 
tained by shifting the blue edge of the bulk of the background 
stars along the reddening curve to the turn-off point (step 1 
in Fig. 1^. Since the majority of the background stars are red- 
dened beyond the selected blue edge and their intrinsic colours 
are significantly bluer than the turn-off point (^MOV), the cor- 
responding extinction offset is a crude lower limit of the real 
value. Although this lower limit is sufficient for most cases, 
it can be further refined (step 2 in Fig. ^ by using knowl- 
edge of the galactic structure and the LF. The Galaxy model 
of Bahcall & Soneira ( |1980| ) was used toge ther w ith the LFs 
in Gould et al. ( |l996| ) {My > 8) and Scalo ( |l986| ) (My < 8) 
to derive the distribution of background stars as a function of 
cloud extinction. Figure |l] shows My at the peak in the distri- 
bution of through-shining stars as a function of Avcioud- The 
observed fields are at a low galactic latitude and suffer from 
high extinction. The background stars are therefore dominated 
by the thin-disk population, which makes the single-component 
disk of Bahcall appropriate to use. Uncertainties inherent in the 
model would mainly arise from systematic errors in the scale 
height vs spectral type relation. Such effects are small enough 
to be dominated by the variation in the LF of the background 
stars. Errors in step 2 that may arise from uncertainties in the 
LF are expected to be small, since the visible background stars 
predominantly consist of dwarfs between F and early K, i.e. 
stars for which the LF is well defined. Steps 1 & 2 are shown 
for the three reference fields in Fig. |[ Note that step 1 is a 
lower limit of the extinction derived independently from any 
assumptions of galactic structure or the LF. The error in step 
2 is estimated to be ^ 0.5 mag. The derived extinction of the 
reference field, Avicf — Ayi + Ayi, is still only a lower limit, 
since we may now use the better estimate Ayrcf instead of Ay\ 
in step 2, and iterate further However, we restrict ourselves 
to two steps and regard the derived value to be a conservative 
lower limit with an error of ^ 0.5 mag in Avrcf ■ 
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Table 1. Observed areas and completeness limits. The photometric system is RI Kron-Cousins (KC) and V Johnson (J) 



field id field centres completeness limits 





RA (2000) 


DEC (2000) 


I 


b 


Ikc 


Rkc 


Vj 


My 


Serpens A 


18:27:29.1 


-3:42:21 


26.8931 


3.6223 


22.1 


23.2 


24.2 


17.2 


Serpens B 


18:28:04.4 


-3:42:21 


26.9610 


3.4921 


22.0 


23.3 


24.4 


17.4 


Serpens C 


18:28:03.3 


-3:49:57 


26.8463 


3.4377 


20.4 


22.3 


23.9 


16.9 


Serpens F 


18:27:52.4 


-3:34:16 


27.0576 


3.5986 


21.3 


22.5 


23.7 


16.7 


Serpens J 


18:28:23.2 


-3:34:18 


27.1164 


3.4847 


21.3 


22.9 


23.4 


16.4 


Serpens K 


18:28:39.4 


-3:45:28 


26.9821 


3.3390 


21.5 


22.9 


24.3 


17.3 


Serpens L 


18:29:13.7 


-3:45:28 


27.0480 


3.2125 


22.2 


23.6 


24.5 


17.5 


p Oph E E 


16:31:56.7 


-24:27:53 


353.8780 


15.9224 


21.5 


22.6 


23.5 


17.6 


p Oph E F 


16:32:49.3 


-24:29:24 


353.9928 


15.7544 


21.4 


22.7 


23.7 


17.8 


p Oph E G 


16:32:48.2 


-24:38:01 


353.8765 


15.6643 


21.5 


22.6 


23.5 


17.6 


p Oph E H 


16:33:22.5 


-24:40:14 


353.9349 


15.5415 


21.5 


22.7 


23.7 


17.8 


p Oph E I 


16:33:24.6 


-24:32:29 


354.0423 


15.6192 


21.5 


22.5 


23.6 


17.7 


p Oph E J 


16:33:23.3 


-24:25:16 


354.1341 


15.7009 


21.4 


22.6 


23.6 


17.7 


p Oph W B 


16:28:00.1 


-24:36:09 


353.1580 


16.5100 


21.4 


22.4 


23.5 


17.6 


p Oph W D 


16:27:42.4 


-24:43:33 


353.0152 


16.4791 


20.3 


21.0 


22.4 


16.5 


p Oph W D' 


16:27:42.3 


-24:28:45 


353.2084 


16.6420 


21.3 


22.5 


23.4 


17.6 


p Oph W E 


16:27:27.7 


-24:36:09 


353.0736 


16.6023 


20.9 


22.6 


23.5 


17.7 


p Oph W G 


16:27:09.1 


-24:27:49 


353.1339 


16.7469 


20.0 


21.5 


22.6 


16.7 


Orion A Fl 


5:41:22.6 


-8:34:35 


212.7511 


-19.4421 


21.5 


23.4 


24.5 


16.5 


Orion A F2 


5:41:22.6 


-8:25:57 


212.6138 


-19.3793 


21.5 


23.4 


24.5 


16.5 



4.1. Orion A 



5. Selecting foreground stars 



The whole area has a high roughly constant extinction. The 
two-step procedure was applied to the entire field, resulting in 
Avref > 6.8. From ^^CO data in Fukui & Mizuno ( |1991| ) an 
independent value of Ay ^ 7 was derived using the ^y^^^CO 



relation in Bally et al. ( 1991 ). 



4.2. Serpens and p Oph 

The fields were divided into overlapping squares of 1' x 1', sep- 
arated by 15". The number of stars within each square was 
summed in 0. 1 mag bins, providing the extinction relative to 
the reference field. 

The results for Serpens and p Oph respectively are ^vrof > 
6.9 and ^vrof > 2.4. The same reference field was used for 
both areas in p Oph, since they are separated by only 1°. 

The extinctions at the centres of the two areas in p Oph 
were estimated from ^^CO column densities in Loren (1989). 
Using the same ^y-^^CO relation as in the Orion case. Ay ^ 
11 was deduced. The present survey gives a lower limit of 
Ay = 7.5. For the reference field, the ^'^CO gave Ay ~ 3, 
whereas this survey gives a lower limit of Ay^^f ~ 2.4. 

For both Orion A and p Oph, the lower limits from the 
present survey are in good agreement with independent values 
derived from ^^CO data. This shows that reliable lower limits 
on the extinction can be obtained directly from the VRI data. 
The Serpens extinction, derived using the same method as for 
Orion and p Oph, is therefore likely to be a reliable lower limit 
as well. 



To really make use of the clouds as screens, crowded regions 
must be avoided. A conservative upper limit on the average 
surface density of stars brighter than the completeness limits 
was set at 4 per arcmin^. The remaining areas are for p Oph 
689 arcmin^ and for Serpens 412 arcmin^, 84 and 77 % of 
the respective original fields. For Orion A the entire field (158 
arcmin^) is sufficiently sparse. 

Figures of the target areas, including selected foreground 
stars and excluded regions are shown in App. A. 

5.7. Colour-colour analysis 

The reddened background stars separate from the foreground 
main sequence at V~R — R—I — 0.9 (MOV), resulting in two 
parallell sequences separated by ^ 0.5 mag mV — R (Fig. 
All stars within 0. 15 mag of the empirical main sequence were 
considered as foreground candidates by colour and had their 
photometric parallaxes measured (see Sect, below). The result- 
ing foreground M-dwarf sample is shown in Figs. ^ - ^together 
with all reddened background stars brighter than the complete- 
ness limit. 

In p Oph a number of bright stars fall inbetween the main 
sequence and the reddened background stars. This is not a re- 
sult of systematic photometric errors, since in the same field 
there are many stars lying perfectly on the reddening sequence. 
These objects are preferably found in p Oph W. Since they are 
absent from the Serpens area at a lower galactic latitude, red- 
dened M giants are not a likely explanation. Most probably they 
are recently formed stars, still embedded in the nebula (see also 
Sect.|andFig.|l0l). 
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Fig. 4. Star brighter than the completeness limit, including the 
final foreground M dwarfs in the Orion A area. Errors are 1 a. 
The dashed line is the main sequence. The solid line is a fit 
to the reddened background stars in p Oph and Serpens. Dou- 
ble ring objects are the finally accepted sample of foreground 
candidates, fitting as such both by colour and photometric dis- 
tance. The two reddest final candidates were also measured in 
J and thereby confirmed to be unreddened objects 



As photometric WD candidates those objects ly ing ap prox- 
imately in the region of the WDs in Bergeron et al. ( |1997| ) were 
considered (Fig. These objects are distinctly separated from 
other stars and form a well-defined isolated group. Note that the 
reddest WDs systematically fall below the empirical sequence. 
This is probably caused by the increasingly different spectra in 
WDs as compared to the standard stars towards cooler temper- 
atures. 



5.2. Photometric parallaxes 

Photometric distances, although suffering from a large spread 
in the Afy vs y — / relation, are used here as an additional tool 
for cleaning out background stars. As a bonus, cloud distances 
get independent consistency checks. Note the lack of stars be- 
tween 400 and 500 pc for Orion A in Fig. ^ clearly indicating 
that d = 480 pc is an overestimate of the distance. 

To check possible foreground extinction, stars within 7° 
from the tar get ar eas were selected from the UBV catalogue 
MermilHod ( 1987| ) via SIMBAD. Eb-v was determined by the 
Q met hod for stars earlier than AO (see e.g. Mihalas & Binney 
( 1981 )). Later stars with known spectral types were also used. 
For all three areas it was concluded that foreground extinction 
is too small to affect the results of this paper. 
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Fig. 5. Star brighter than the completeness limit, including the 
final foreground M dwarfs in the Serpens area. Symbols are as 
in Fig. ^. Single rings are foreground candidates by colour, but 
not by distance 
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Fig. 6. Star brighter than the completeness limit, including the 
final foreground M dwarfs in the p Oph area. Symbols are as in 
Fig.|&| 



5.2.1. M dwarfs 

Each object was assigned a mean, min and max distance. The 
mean distance is derived from the composed My V — I 
relation in RG97. Min and max values come from their cor- 
responding 1 a spread . A histogram of the mean distances to 
the colour-colour foreground candidates marked in Figs. 0- |]is 
shown in Fig. ^. Stars consistent with being foreground objects 
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Fig. 7. Photometric WD candidates. Symbols are as in Figs. ^- 
^. New sy mbols : The shaded area outlines the WDs of Berg- 
eron et al. (|1997|). Errorbars are 1 cr. J-1 13 is a WD from JDH94 
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Fig. 8. Photometric mean distances to the photometric fore- 
ground M-dwarf candidates. The solid bars show the adopted 
cloud distances 



are distinctly sorted out. The final sample of probable fore- 
ground M dwarfs consists of 21 stars (Table ^). 

5.2.2. White dwarfs 

Photometric distances were estimated via a linear least-square 



fit of A/v vs F - / of the WDs in Bergeron et al. ( |1997[ ) with 
parallax errors less than 30 % : 
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Fig. 9. Low-resolution spectra of two of the Serpens WDs. The 
dotted lines show Hel lines, dased lines are Hell and solid lines 
H. Note the sky lines at AA5890, 6867 (Na D, and atmospheric 
B band) 
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This functional form provides only a rough estimate of the 
WD absolute magnitude. One should keep in mind that there 
are WDs like ES0439-26 (Ruiz et al. |1995| ) that do not follow 
this relation in the least. 

It is possible, although unlikely, that these WD candidates 
are reddened background stars. Since the background extinc- 
tion at the position of the WD candidates in all cases is esti- 
mated to exceed 6.8 mag in V (Table a through-shining star 
of the same colour as our reddest WD candidate {V ~ I ^ 1.5) 
would have an unreddened colour {V — I)o < —1.2. There 
are no such stars. Suppose that they are indeed early-type stars 
shining through low-extinction windows. A BIV star behind 
OrionA would then have to be at a distance of about 0.5 Mpc 
to fit the observed magnitudes, i.e 160 kpc above the galactic 
plane, which is unrealistic. The simplest answer is that these 
objects are indeed foreground WDs. 

5.2.3. Spectroscopic confirmation of white dwarfs 

In Fig. ^ spectra taken at the NOT of the two brightest WDs are 
compared with two spectrophotometric standard stars. Hz21 is 
a He-rich WD of type DO, while G93-48 is of type DA. The 
serpwd3 spectrum is classified as typ e DB, since it only shows 
Hel lines (Dreizler & Werner 1996). The second WD is ser- 



pwdl, showing a rather noisy spectrum due to its faint mag- 
nitude. However, the absence of prominent spectral features 
in combination with its magnitude and colour makes it a very 
Ukely WD. 

6. Constraints on foreground candidates from other 
surveys 

The Orion A area (LI 64 1) has been extensively su rveye d in 
the infrared (Strom et al. (|l989l |l993|); Ali & Depoy |l995|). In 



L. Festin: The luminosity function of white dwarfs and M dwarfs 



7 



Table 2. Final foreground candidates. Photometric errors are internal and 1 a. The photometric system is Rl Kron-Cousins (KC) 
and V Johnson (J). Upper and lower distance limits are also 1 cr (cr originating mainly from cosmic spread in the My - V ~ I 
relation. Ay is the lower limit of the background extinction at each star's position. The last column is the weighted count each 
star gets when reduced to the solar neighbourhood 



id 




RA (2000) 


DEC (2000) 


Ikc 


Rkc — Ikc 


Vj - Ikc 


Mv 


dm 


d+ 


d- 


Ay 


w.ct. 


omd 1 




5 


41 


22.7 


-8:22:21 


15.77 0.01 


1.50 0.01 


2.64 0.01 


11.6 


230 


195 


270 


6.8 


1.32 


omd 2 




J 


A 1 

41 


JO.O 


-8:30:18 


15.95 0.01 


1.26 0.01 


2.29 0.01 


10.6 


346 


295 


407 


6.8 


1.32 


omd 3 




5 


41 


14.0 


-O.JO, ly 


1 ^ QQ n n 1 


1 n ni 


Z.JO U.Ul 


1 n s 

lU.O 




Z /O 


IS 1 

JOl 


O.o 


1 ^9 
L.JZ. 


omd 4 




5 


41 


17.5 


-8:29:44 


16.10 0.01 


1.50 0.01 


2.65 0.01 


11.7 


267 


228 


314 


6.8 


1.32 


omd 5 




5 


41 


04.5 


-8:38:23 


17.26 0.01 


1.53 0.01 


2.74 0.01 


12.2 


381 


308 


471 


6.8 


1.32 


omd 6 




5 


41 


14.8 


-8:25:45 


17.34 0.01 


1.58 0.01 


2.89 0.01 


13.0 


290 


235 


358 


6.8 


1.32 


omd 7 




5 


41 


19.8 


-8:30:48 


17.87 0.01 


1.62 0.01 


2.89 0.01 


13.0 


370 


300 


458 


6.8 


1.32 


omd 8 




5 


41 


31.8 


-8:22:39 


18.09 0.01 


1.83 0.01 


3.27 0.01 


14.3 


262 


232 


295 


6.8 


1.32 


omd 9 




5 


41 


23.7 


-8:35:42 


18.35 0.01 


1.81 0.01 


3.21 0.01 


14.2 


300 


266 


338 


6.8 


1.32 


omd 10 




5 


41 


21.5 


-8:23:09 


19.63 0.02 


2.00 0.02 


3.72 0.06 


15.5 


382 


339 


430 


6.8 


1.32 


omdl 1 




5 


41 


18.9 


-8:27:39 


20.12 0.01 


2.33 0.04 


4.17 0.11 


16.6 


358 


317 


403 


6.8 


1.32 


serpmdl 


18 


27 


10.8 


-3:43:59 


16.66 0.01 


1.64 0.01 


2.94 0.01 


13.2 


186 


151 


230 


7.6 


1.03 


serpmd2 


18 


27 


57.2 


J.J i . J / 






Z..Z. 1 \J.\J 1 


lyJ.J 


241 


70'^ 


9X^ 

^OJ 


1 1 ^ 


1 0^ 


serpmd3 


18 


27 


53.8 


-3:36:53 


16.10 0.01 


1.68 0.01 


3.02 0.01 


13.6 


125 


110 


140 


8.9 


1.03 


serpmd4 


18 


28 


32.8 


-3:34:53 


16.81 0.01 


1.61 0.01 


2.91 0.01 


13.1 


212 


171 


262 


9.9 


1.03 


serpmdS 


18 


29 


20.3 


-3:44:15 


17.29 0.01 


1.66 0.01 


3.03 0.01 


13.7 


212 


189 


240 


12.1 


1.03 


ophmdl 


16:33 


02.6 


-24:30:14 


14.60 0.01 


1.54 0.01 


2.84 0.01 


12.7 


88 


71 


108 


4.8 


1.09 


ophmd2 


16:32 


42.8 


-24:26:55 


15.93 0.01 


1.67 0.01 


3.08 0.01 


13.8 


109 


97 


123 


6.6 


1.09 


ophmd3 


16:27 


30.4 


-24:32:34 


14.43 0.01 


1.38 0.01 


2.48 0.01 


11.1 


143 


122 


168 


7.8 


1.09 


ophmd4 


16:27 


36.5 


-24:28:33 


13.46 0.01 


1.42 0.01 


2.47 0.01 


11.1 


92 


78 


108 


7.8 


1.09 


ophmdS 


16:27 


26.6 


-24:25:54 


15.99 0.01 


2.09 0.01 


3.93 0.01 


16.0 


60 


53 


68 


7.8 


1.09 














Ikc 


Rkc 


Vj 














owd 1 




5:41:37.7 


-8:36:09 


22.00 0.14 


22.78 0.06 


23.25 0.04 


15.8 


313 


269 


365 


6.8 


1.32 


owd 2 




5:41:04.0 


-8:23:09 


21.24 0.03 


21.89 0.03 


22.38 0.02 


15.5 


238 


207 


273 


6.8 


1.32 


owd 3 




5:41:22.6 


-8:22:50 


22.83 0.10 


23.55 0.15 


24.08 0.12 


15.8 


459 


393 


537 


6.8 


1.32 


serpwdl 


18:28:05.6 


-3:40:02 


20.38 0.02 


20.79 0.01 


21.14 0.01 


14.6 


208 


181 


239 


9.6 


1.03 


serpwd2 


18:27:45.0 


-3:46:01 


20.51 0.02 


20.98 0.02 


21.47 0.01 


15.0 


192 


167 


221 


9.9 


1.03 


serpwd3 


18:27:41.6 


-3:40:38 


17.88 0.01 


17.90 0.01 


17.87 0.01 


12.6 


111 


96 


127 


12.6 


1.03 


serpwd4 


18:28:26.3 


-3:43:52 


21.14 0.05 


21.96 0.04 


22.61 0.02 


16.3 


181 


157 


208 


11.6 


1.03 



" JciT = 17.77 (0.06) 
* JciT = 18.08 (0.07) 



Strom et al. (1989) there is one overlapping object, however 
not a foreground one. The other two surveys are deeper but 
unfortunately do not overlap. 

In p Oph, the western area was entirely covered in JHK 
in a survey by Barsony et al. (1997). There are 30 overlap- 
ping sources brighter than our completeness limits, including 
all three foreground candidates in this area (ophmd3-5 in Ta- 
ble In Fig. |l^ the appearance of the overlapping sample 
in VRI is compared to UK. Ophmd3 and 4 are confirmed as 
foreground by the additional JHK data. For ophmdS, the three 
independent JHK measurements are not consistent, implying 
either variability or a problem in the photometry. The fact that 
ophmdS does not keep its position relative to the bulk of the 
background stars between VRI and UK supports it being a 
foreground star of spectral type about M6V. 



The p Oph E area was covered in VRI by Jarrett (JDH94 
-I- thesis). Both our candidates in this area (ophmdl-2) have 
magnitudes consistent with theirs, and are considered as pos- 



sible foreground stars in Jarrett (1992). Ophmd2, though, was 



rejected from the foreground candidates in JDH94. Both stars 
are slightly redder than the main sequence mV — R (Fig. 
but are still regarded as foreground stars. See Sect. |]l| for a 
thorough comparison to Jarrett's results. 

In our Serpens field, no complementary photometric data 
could be found. 



7. The luminosity function 



A simplified version of the Vmax method (Felten 1976) 



outlined below, was used to derive the LF. In a photometric 
survey each star has a maximum distance, dmax, at which it 
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□ ''<] 



area 


Vmax 


Vmax,gcn 


Orion A 


285 


216 


Serpens 


193 


188 


pOph 


66 


60 


total 


544 


464 



M6V 
'M5V 



M3V 



3 

l-J 



Table 4. The foreground M-dwarf sample binned into 1-mag 
bins in My. The 2nd last bin is complete in 70 % of the volume 
and the last bin in 26 %. No correction for incompleteness has 
been applied 







bin 




systems 




binary corrected 






My 


iV 


A^gcn 


cum 


iVgcn 


cum 






10-11 


3 


3.67 


3.67 


4.49 


4.49 






11-12 


4 


4.82 


8.49 


6.19 


10.68 






12-13 


2 


2.41 


10.90 


4.46 


15.14 






13-14 


7 


7.85 


18.75 


9.30 


24.45 






14-15 


2 


2.64 


21.39 


4.20 


28.65 






15-16 


1 


1.32 


22.71 


2.62 


31.27 




.---'''m5V • 


16-17 


2 


2.41 


25.12 


3.06 


34.33 


_.---Gr---""' 
M3V 




17-18 








25.12 


0.80 


35.13 







1.4 



1.8 



2.2 



2.6 



R-l 



Fig. 10. Overlapping near-IR sources. The two earliest fore- 
ground candidates are represented in the upper panel by mean 
values of the JK magnitudes given i n Bar sony et al. (1997), 
Strom et al. ( 1995 ) and Greene et al. ( 1992 ). Individual values 
were used for the latest candidate, triangle ( Barsony), penta- 
gram (Strom) and box (Greene). Dashed lines are the corre- 
sponding main sequences and the solid lines show the empiri- 
cal reddening sequences 



still can be detected. The corresponding maximum volume, 
Knaxi is defined by the survey area and dmax- Each star's 
contribution to the LF is then V^^^^. In order to compensate 
for galactic density gradients and reduce the derived LF to the 
solar neighbourhood, Vmax must be replaced by Knax.gcn, a 
generalized volume (Tinney 1993): 



- ^2^-rsin\b\/h^^ 



where il denotes the angular area of the field, h the disk 
scale height and b the galactic latitude. The distances involved 
here are small enough to neglect the radial density gradient 
in the Galaxy. Both the M-dwarf and WD distribution was 
approximated by the old disk population, with a scale height 
h ~ 350 pc (Bahcall & Soneira 1980). Deriving V^max is a 
trivial operation in a dark cloud survey, since for all stars 
brighter than the completeness limit (in absolute magnitude) 
rfmax = c^cioud- This is an important point of this paper, the 



survey volume not being derived through a colour-magnitude 
relation having a large cosmic spread. Table || gives the 
physical and generalized values of the foreground volumes. 
The last column in Table ^ shows the weighted count for 
each foreground star deduced as Knax/Vmax.gcn- Errors in the 
distances to the clouds are expected to be on the order of 10 %, 
implying volume uncertainies of 30%. This uncertinty has 
not been included in the further analysis. 



7.1. M dwarfs, main sequence 

The majority of foreground stars are early-type M dwarfs, 
2.5 < V — I < 3.2, reflecting the well-known peak in the 
LF. The low number of M dwarfs (21) implied a rather coarse 
binning in the LF. Table ^ gives the LF in 1-mag bins, while 
in Fig. |ri]2-mag bins were chosen. The agreement with previ- 
ously determined LFs is satisfactory. 

For stars brighter than My — 10 (~ Ml V) the number that 
is lost due to saturation starts to become significant. These stars 
are therefore not considered further. 

The large spread in the faint bins in Fig. |ll| illustrates the 
present uncertainty in the stell ar LF. The Poissonian 1 <t uncer- 
tainties in Jahreiss & Wielen ( |l997| ) and RG97 are about 30 % 
in the last two bins. 



7.1.1. Lost companions 

RG97 found that the number of lost unresolved companions 
approximately is compensated for by distant stars spread into a 
magnitude limited sample via Malmquist bias. This will not be 
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60 
40 
20 



^ 24 
> 



16 



cumulative LF 



O 

-0- 
it 



normal LF 



* 
O 



10 



11 



12 



13 



14 
M,, 



16 



17 



18 



2.4 
2.2 
2 



1.8 

> 



1.6 
1.4 



1.2 





J21f 












U221 




J74 



J43 



1.4 



1.6 



J19,-- J220 



1.8 2 
R-l 



2.2 



2.4 



Fig. 12. A comparison with JDH94. The dashed Hne marks the 
main sequence and the dots outHne the reddening sequence. 
The + signs mark the overlapping objects with 1 a errorbars. 
The JDH94 counterparts are marked with Jx, x indicating the id 
in Table 3 in JDH94. The dotted lines show the deduced colour 
differences. For J74 there are two independent measurements 
in our survey 



Fig. 11. The M-dwarf luminosity function. The dashed line is 
the system LF and the solid line the single-star LF derived af- 
ter binary correction. Diamonds and stars show the LF of Kirk- 
patrick et al. (1994) and Gould et al. (1996) respectively. As- 
terisks and circles outline the LF of the nearby stars in RG97 



and Jahreiss & Wielen (1997) respectively. The errorbars show 
Poisson errors of the generalized counts before binary correc- 
tion. Errors in the LF arising from uncertain cloud distances 
have not been included (expected to be ~ 30%). The incom- 
pleteness in the last bin (see Table ^ has not been corrected 
for 



the case in this survey, since the foreground sample is strictly 
volume limited and there is a dark cloud preventing distant stars 
from interfering. The number of lost companions was estimated 
from the binary distribution in the nearby stars in RG97. After 
applying the correction for galactic structure 10 companions 
were expected to have been lost. These were distributed as the 
lost binary companions in Fig. 12 of RG97 and added to our 
system LF. The resulting companion-corrected LF is shown in 
Fig. [nl and Table |. 

7.2. White dwarfs 

The sample of WDs is too small to derive any details in the WD 
LF. Seven WDs were found, 0.015 (±G.006)pc~^ (Poissonian 
errors), which is a significantly higher number than expected 
from LDM88 (0.0 03 pc ^), but in better agreement with RT95 
and Oswalt et al. ( |1996| ), who both give ~ 0.008 pc"^. These 
numbers will be dicussed in more detail in Sect. 8.2. 



8. Dicussion, comparison to other works 

For the M-dwarf part, the agreement with other photometric 
surveys is satisfactory (Fig. |ll]). It is clear that any contamina- 
tion of the M-dwarf sample from background stars and/or stars 
associated with the clouds must be small and that the opaque- 
nebula method works reliably also in a purely photometric ver- 
sion, provided that the photometry is sufficiently accurate. 

The main advantage in comparison to other methods is 
that the survey volume does not rely on photometric paral- 
laxes. Even though the individual absolute magnitudes do suf- 
fer from cosmic spread in the colour-magnitude relation, each 
star's contribution to the space density does not. 

If the distance to the cloud is not known a priori, it can 
either be estimated from the photometric distances of the fore- 
ground objects or by their integrated number in a well-known 
part of the LF. In this context it would be useful to obtain more 
data for foreground stars brightwards of the peak in the LF at 

My ^ 12. 

The derived samples are volume limited and unaffected 
by dynamical bias. Finally, the dark clouds effectively screen 
background galaxies and giant stars from the foreground pop- 
ulation. 

Among the drawbacks are that the selected foreground ob- 
jects in general will be too distant for binaries to be resolved 
and therefore the derived LF must be corrected for by an a pri- 
ori known binary distribution function. The dark clouds are not 
evenly spread across the sky, which means that the available 
areas are restricted. On the other hand, they do provide a tool 
for probing areas near the Milky Way plane. 
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8.1. The Jarrett paper 

We compare our result with that of JDH94. The generalized 
volume sampled in that paper is 52 pc^ (p Oph) +101 pc'' (Tau- 
rus), totally 153 pc^. A striking difference as compared to our 
survey is their excess of late foreground candidates, seen also 
in the Barnard 5 survey (Jarrett 1995 ). The comparison here 
was made in p Oph E, where part of our surveys overlap. Five 
of their six foreground candidates that were cross identified lie 
perfectly on the reddening sequence in our data and are defi- 
nitely not foreground objects (Fig. 12). Since our data for the 
bright J43 are consistent, the photometric calibration is likely 
to be correct. A possible explanation of the colour discrepancy 
is that these sources were spread into the JDH94 foreground 
sample as a result of their lower signal-to-noise ratio. Our data 
reach ^ 2 magnitudes deeper than JDH94. It may also be that 
the JDH94 photometry is subject to systematically discrepant 
corrections for atmospheric correction, since their data were 
acquired at high air masses (^ 2, Jarrett, priv comm.). The ten- 
dency that very faint stars, barely measurable, scatter into the 
M-dwarf zone is seen also in our data, and was the main reason 
for including objects only to the completeness limit and not all 
the way down to the detection limit of the survey, about 1 mag 
fainter. 



8.2. White dwarfs 

A local number density of 0.015 (±0.006) WDspc^'^ was de- 
rived from the 7 WDs discovered in this survey. This is a fac- 
tor 5 more than predicted from LDM88. RT95 found 0.008 
(±0.003) WDs pc"'^ in a proper motion survey, consistent with 
our result, but as in our case the numbers are small. An interest- 



The results of the MACHO micro lensing experiment (Al- 
cock et al. 1997 1 have given a mean mass of the lensing objects 



ing result is that Ruiz et al. ( |1993| ) found a factor 2 more stars 
with II > 0.5" yr~^ than present in LHS (Luyten 1976| ) for 
overlapping areas. The LDM88 sample was drawn from LHS 
and would also be affected by this proposed incompleteness, a 
fact that may explain the discrepancy in the LDM88 WD den- 
sity as compared to more recent results. 



Oswalt et al. ( 1996 ) used yet another approach and arrived 
at 0.0053 j:o"!^gg7 WDspc-3 for WDs in wide binaries. They 
added the single-WD LF in LDM88 (0.0023 WDspc-^) and 
concluded that the total space density of WDs in the solar 
neighbourhood is 0.0076to"ooo7 WDspc"^. Our results and 
RT95 show that the LDM88 sample is likely to suffer from 
a rather large incompleteness, suggesting that adding 0.0023 
WDspc^'^ to the binary part results in a substantial underesti- 
mation of the total WD space density. 

Our sample of WDs includes neither wide nor narrow bina- 
ries, except possible double WDs. If we add our WD density to 
the Oswalt binaries the result is 0.02 (±0.007) WDs pc~3, cor- 
responding to 0.013 (±0.005) Mq pc^'^ if we assume a mean 
mass of 0.65 Mq per WD. The WD fraction of the dynami- 
cal mass would then be 17 (±7)%. Object s like ESO 439-20 
(RT95) and WD 0346+246 (Hambly et al. 
the WD density may be even larger. 



of ^0.5 Mq, although strongly model dependent. The MA- 
CHO mass density in the solar neighbourhood is on the order 
of 0.005 Mq pc"'^, which in fact is consistent with the pro- 
posed increase of the WD space density. However, more defi- 
nite results will have to await a larger volume limited sample 
of WDs. 



9. Conclusions and future prospects 

The space density of M dwarfs and WDs has been studied as 
volume limited samples in front of dark high extinction nebu- 
lae, a method that efficiently reduces the impact of Malmquist 
bias, since the survey is volume limited and the LF not de- 
pendening on individual photometric distance estimates. Areas 
towards the molecular cloud complexes p Oph, Serpens and 
Orion A were observed. The survey covers a volume corre- 
sponding to 464 pc'^ in the solar neighbourhood, complete to 
A/v ^ 16.5. In total 21 M dwarfs and 7 WDs were identified 
as foreground objects. 

Although details are hard to resolve due to the limited sam- 
ple, our M-dwarf LF is consistent with previous photometric 
results, and after binary correction also with the nearby stars. 
Our survey does not show any signs of the reported upturn in 
the fain test magnit ude bins in th e prev ious dark cloud surveys 
(Jan-ett |l992[|l995| ; Jarrett et al. |l994| ), and we conclude that 
their excess of faint M dwarfs is probably caused by their larger 
photometric errors. 

The unexpected appearance of as many as 7 WDs com- 
pared to the 1.4 estimated from the WD space density in 
LDM88 is perhaps the most interesting result of this paper. 
The derived space density is 0.015 (±0.006) WDspc^'^. By 
combining our sample with the binary WDs in Oswalt et al. 
( |1996t ) the final estimate is 0.02 (±0.007) WDs pc-^, implying 
a mass density of 0.013 (±0.005) Mqpc^'^. This means that 
the WDs may contribute by as much as 15-20% to the local 
dynamical mass. No doubt this subject is not yet settled and 
deserves further investigation. 



1997) suggest that 



The method of opaque nebulae has indeed proved to be 
efficient in identifying foreground M dwarfs and WDs. Con- 
cerning the latter, any nebula that provides a colour excess 
Ev-i > 2 (Ay > 5) ensures that also the bluest main se- 
quence stars get a reddened colour beyond the reddest known 
WDs (V — I ^ 1.5). Already for Orion A the available area is 
■~ 50 times larger than our field, enclosing a volume that would 
host 150 WDs and ^^550 M dwarfs, numbers that certainly 
are sufficient for a more detailed analysis of the corresponding 
LFs. The increased availability of large field CCDs and CCD 
mosaics provides the ideal tools for this kind of survey. We 
have also shown that the distance to the clouds and the cloud 
extinction may be reliably derived from the VRI data alone, 
an important point if a future large-scale survey of this kind is 
to be carried out. 
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Fig.A.2. Serpens 4° south of the Serpens cloud core). Coordinates are equinox 2000. Symbols as in Fig. 4.1. Filled areas 
were rejected from the survey due to high stellar density (Sect. H) 
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Fig. A.3. p Oph W (left, L1688) and p Oph E (right, L1689). Coordinates are equinox 2000. Symbols as in Figs. A. 1 - \.2 
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Fig. A.l. The Orion A area (southern part of L1641). Coordi- 
nates are equinox 2000. Small dots mark all objects observed 
in /. Rings are the foreground M dwarfs. Pentagrams are white 
dwarfs 



